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ABSTRACT 


A rapid analytical technique is used to deter- 
mine simultaneously the possible configurations of an 
elliptical parking orbit at Mars. Both the geometry 
of planetary approach and departure and the effects of 
planetary oblateness are used to determine the char- 
acteristics of those orbits which shift into proper 
alinement for departure and require no discrete pro- 
pulsive maneuvers. The technique makes use of the 
separable nature of the equations that express secular 
rotation of the orbit plane and major axis. Prelimi- 
nary results indicate that the orbital inclinations fall 
into narrow bands which are widely distributed over 
the interval of possible inclinations. 


USE OF PLANETARY OBLATENESS FOR 


PARKING-ORBIT ALINEMENT 

By Joseph R. Thibodeau HI 
Manned Spacecraft Center 


SUMMARY 


A technique relying upon the effects of gravitational harmonics for parking-orbit 
alinement is analyzed. The basic approach is to choose orbital eccentricity and incli- 
nation so that the resulting nodal and apsidal motions will shift the original orbit into 
proper alinement for departure on the intended date. Thus, the need for propulsive 
maneuvers to change the orbital inclination, line of apsides, or orbital eccentricity are 
reduced or eliminated; and periapsis impulses can be used for both the capture and the 
escape maneuvers. Numerical results are presented for a long- stay-time Mars mis- 
sion in 1977. These results illustrate the proposed technique, but do not include com- 
parisons with alternative methods. 


INTRODUCTION 


During planetary capture and landing missions, long-duration parking orbits may 
be established about the target planet. These stopover missions often require substan- 
tial, and sometimes costly, readjustment of the parking orbit to provide for efficient 
planetary escape and return to Earth. Rotation of the parking orbit becomes necessary 
to compensate for two dynamic effects: (1) the orientation and motion of the hyperbolic 

departure asymptote and (2) the motion of the parking orbit resulting from the aspher- 
icity of the planetary gravitational field. 

Efficient propulsive techniques are available to accommodate the two dynamic 
effects (ref. 1). These techniques require auxiliary maneuvers, which may be either 
discrete or combined with the planetary capture and escape maneuvers, to accomplish 
orbital readjustment. 

Nominal missions which require no auxiliary maneuvers can be designed by rely- 
ing on the effects of the gravitational harmonics. This technique was originally pro- 
posed in references 2 and 3. An analytical approach to the technique will be discussed. 
This technique accommodates both the geometry of approach and departure and the 
effects of planetary oblateness. It also permits a variety of elliptical parking orbits 
to be simultaneously found. 



The technique is intended for application in preliminary trajectory design and 
early feasibility studies of orbital missions. In these instances, analytical techniques 
can be used advantageously to evaluate many possible configurations of the parking 
orbit. The method is particularly useful for the study of eccentric parking orbits which 
match so-called " double -Hohmann" transfers, or minimum- energy members of the 
conjunction class of Mars missions. Conjunction- class missions generally have total 
trip times of 950 to 1000 days with stay times ranging from 300 to 500 days. For these 
missions, there is sufficient time for perturbations to rotate high-eccentricity parking 
orbits through large angles. This method is somewhat less successful for the study of 
parking-orbit configurations required for the opposition-class missions in which stay 
times on the order of 40 days are investigated. For these short stay times, fewer 
parking-orbit configurations are available, and the parking orbits generally have low 
eccentricities. 

Only secular rotation of the orbital plane and major axis due to planetary oblate- 
ness is considered. Orbital realinement is produced by rotation of the orbit line of 
nodes and the line of apsides, and these rotations are primarily the result of planetary 
oblateness. In the presence of a third body, there are long-period rotation rates caused 
by coupling of the oblateness and third-body perturbations (described in refs. 4 and 5). 
In addition, resonances can lead to feedback between the orbit and the perturbing func- 
tion (ref. 6). 

In general, inclusion of a third body and the higher order harmonics of the plane- 
tary gravitational field requires numerical integration. The technique developed in this 
report can be used as a tool for preliminary mission planning. This technique permits 
rapid determination of all possible parking-orbit configurations. Orbits which best 
meet the required ground rules or boundary conditions can then be selected for further 
study using precision techniques. 

The author wishes to thank Victor R. Bond for the use of the interplanetary 
trajectory-analysis program for the example mission. 
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mean motion of the parking orbit 
planetary equatorial radius 
radius of periapsis 

parking-orbit regression time 

stay time in parking orbit 

rotation time of orbital node 

rotation time of periapsis vector 

hyperbolic excess velocity 

unit vector of hyperbolic excess velocity 

coordinates in the planetocentric system 

planetocentric longitude 

impulsive velocity 

planetocentric declination 

auxiliary angle defined by equation (4) 

gravitational parameter 

auxiliary variable = sin '*‘(tan 6 /tan i) 
ascending node 

unit vector in the direction of the ascending node 
normalized rotation rate 

nodal regression rate 

argument of periapsis 

periapsis precession rate 
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Subscripts: 

A 

D 

h 

P 

S2. 


n. 


approach trajectory 
departure trajectory 
hyperbola 
periapsis 

ascending node closer to the vector, defined by equation (1) and 
figure 2 

ascending node farther from the vector, defined by equation (1) 
and figure 2 


ANALYSIS 


Technique for Orbital Alinement 

The sequence of maneuvers for the orbital-alinement technique is illustrated in 
figure 1 and is outlined as follows: 

1. Orbital insertion occurs at periapsis of the arrival hyperbola. At the time of 
periapsis passage, the following conditions exist: 

a. The parking orbit and the arrival hyperbola are coplanar. 

b. The periapsis position vectors of the parking orbit and arrival hyperbola 
are identical. 

2. During the orbital coast phase, orbital perturbations continually rotate the 
orbital plane and major axis, thereby shifting the parking orbit into alinement with the 
nominal departure asymptote. 

3. At the nominal time of departure, orbital launch occurs at periapsis of the 
parking orbit, and the following conditions exist: 

a. The parking orbit and the departure hyperbola are coplanar. 

b. The periapsis position vectors of the parking orbit and departure hyper- 
bola are identical. 

The technique implicitly defines two constraints or boundary conditions; conti- 
nuity of the periapsis position and of the orbital inclination must be enforced at the 
match points of the parking orbit and the arrival and departure hyperbolas. For this 
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analysis, radius of periapsis and orbital 
inclination can be regarded not only as 
continuous but constant, since planetary 
oblateness does not contribute secular 
rates of change to either parameter. 

The boundary conditions for the 
maneuver sequence are defined by four 
quantities which are the necessary prob- 
lem input data: 

1. Radius of periapsis r^ 

2. Orbital stay time T 

s 

3. Hyperbolic excess velocity of 
the approach trajectory ^ 

4. Hyperbolic excess velocity of 
the departure trajectory ^ 



Figure 1. - Geometry of the parking orbit at 
the time of arrival and departure. 


These data are calculated by an independent trajectory-analysis program. The inter- 
planetary trajectory is therefore assumed to be known, and only those operations in- 
side the planetary sphere of influence are considered. 


Derivation 

The following parameters are used to define the parking orbit: 

1. Inclination i 

2. Eccentricity e 

3. Radius of periapsis r p 

4. Longitude of the ascending node a ^ 

5. Argument of periapsis o> 

r 

The analysis shows how these parameters are found for regressing parking orbits 
which shift into proper alinement for departure. The orbital parameters are refer- 
enced to the time of periapsis passage on the arrival hyperbola. 

The approach relies on the separable nature of the equations which express secu- 
lar rotations of the orbit plane and major axis. These secular rotation rates are func- 
tions of inclination and eccentricity. The technique allows the inclination function to be 
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separated from the eccentricity function, thereby permitting calculation of the required 
orbit inclination independently of eccentricity. 

The determination of the orbital parameters thus involves two distinct phases. 
The first phase is a search for the orbital inclination, ascending node, and periapsis 
vector. The orbital inclination is found by choosing it as the independent search vari- 
able. By sweeping through the range of possible inclinations, a parametric search 
table is generated, and the proper orbital parameters are then found by interpolation. 
The second phase is a search for the orbital eccentricity. The eccentricity is found 
by using the Newton- Raphson method. 

Planetocentric coordinate system . - It is assumed that the hyperbolic excess 
velocities are referenced to an inertial planetocentric coordinate system. Since inter- 
planetary trajectories are generally calculated in a heliocentric reference, these ve- 
locities must be transformed into this coordinate system. The positive Z-axis of the 
planetocentric coordinate system points north along the planetary axis of rotation. The 
XY-plane is defined by the planetary equatorial plane. The positive X-axis is defined 
by the intersection of the planetary orbit and equatorial planes and corresponds to the 
descending node of the planetary orbit on the equatorial plane. The positive Y-axis is 
90° east of X and completes a right-handed system. 

Further discussion of this coordinate system and a presentation of Mars trans- 
formation equations are provided in reference 7. The position of the North Pole and 
the physical ephemeris of Mars are also discussed in reference 8. 

Basic equations . - Assuming that the respective hyperbolic excess velocities and 
the radius of periapsis are specified, the following equations define the planetocentric 
longitude of the ascending node and the argument of periapsis of both the planetary cap- 
ture and the planetary escape hyperbolas. These data will be used later to determine 
the inclination of the parking orbit. 


The node and periapsis position vectors are double-valued functions of inclination. 
This can be seen in figure 2, which shows 
that two orbital planes contain the 


vectors for all inclinations. 

The longitude of the ascending 
node is determined from 


a n = a - a 
”l 


o 0 = a + <7 + IT 
°2 


( 1 ) 


where o = sin *(tan 6/tan i), and a and 
6 are the planetocentric longitude and 



Figure 2. - Geometry of trajectory planes. 
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declination of the vector. As shown in figure 2, a ^ is the longitude of the 

ascending node closer to the vector. Likewise, is the longitude of the node 

2 

farther from the V vector. The declination of the \T vector defines the minimum 

OO 

orbital inclination, since a is not defined for i less than 6 . 

The semimajor axis of the approach hyperbola a^ is determined from the vis 
viva equation 


a h = 



( 2 ) 


When the radius of periapsis is specified, the eccentricity of the approach hyper- 
bola is 


e 


h 



(3) 


The angle between the arrival V M vector and the periapsis of the approach hyperbola 
can now be found by 


9 = cos 



(4) 


where 


e^ is the eccentricity of the approach hyperbola. 

The argument of periapsis co for the approach hyperbola is 

lr 




P 


cos -1 (n • v^) 


- 6 


(5) 


At the time of periapsis passage, the argument of periapsis of the parking orbit is the 
same as that for the approach hyperbola. 
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Equations (1) to (5) are used to find the node and periapsis vectors at the time of 
injection onto the departure hyperbola, with the following exception: 


6 = cos 



( 6 ) 


where e^ ^ is the eccentricity of the departure hyperbola. 

Determination of orbital inclination . - The longitude of the ascending node 

and the argument of periapsis are now assumed to have been calculated for the 

node closest to the V M vector for both the arrival and the departure hyperbolas, and 

the parking orbit shown in figure 1 is considered once again. During the planetary stay 
time, the ascending node of the approach hyperbola must move to the ascending node of 
the departure hyperbola, and r . must move to r 

p, A p, Lf 

The angular distance the parking-orbit node must travel is 


= « 


a. 


C2,D n,A 


= f(V_ 


,A- V - D'V‘) 


(7) 


where the subscripts A and D indicate the arrival and departure hyperbolas, 
respectively. Likewise, the angular distance the periapsis vector must travel is 


p - "p, D - "p, A - K V “, A’ V «, D’ r p- *) 


Aw = 


(8) 


The speed of rotation of the node £2 and periapsis w because of secular vari- 

s s 


ation Jg is given in reference 9 as 


-3nJ 2 lT 

0 = s cos i 

S rt / rt \ ^ 


2a 2 (l - e 2 ) 


(Jj 

s 


-3nJ 2 R 2 
2a 2 (l - a 2 )' 


(5 . 2. 
2 SU1 1 



(9) 


( 10 ) 
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The time for the parking-orbit node to traverse Aa^ is 


Aa 




S2 




(ID 


and the time for the periapsis to traverse Aw^ is 


Aco 

t = I 

W d) 


( 12 ) 


If the stay time in the parking orbit is designated by T , the problem is to obtain 

s 

the inclination and eccentricity of the parking orbit such that 


T s “ " t U ) 


(13) 


This condition is obtained in two steps: (1) A parking-orbit inclination that will make 

t fi equal to t^ is found, and (2) the parking-orbit eccentricity is adjusted to make t^ 

and t^ equal to the planetary stay time T . 


It should be noted that a closed-form solution is not evident. If equations (11) and 
(12) are equated, a solution for the inclination would then be found in terms of Aa fi 

and Ao> . But Aa„ and Aw are not known, since they are likewise functions of 

P 06 p 

orbital inclination, as is shown by equations (7) and (8). 


A simple interpolation scheme lends itself to finding the orbital inclination for 
which t^ equals t^. The basic equations are used to calculate t^ and t^ at 10° 

increments through the interval of possible inclinations. The ratio can be tab- 

ulated with inclination, and the equality of t^ and t is indicated by a ratio of 1. 

Thus the table can be scanned to find all the time ratios near 1, and the possible 
parking-orbit inclinations can be found by interpolation. If necessary, the previously 
discussed procedure can be repeated, using a smaller inclination increment to make 
t^ and t equal within some specified tolerance. 


In practice, not one, but four ratio-versus- inclination tables are simultaneously 
generated. As shown in figure 2, there are two trajectory planes containing the 
approach asymptote for every inclination. Thus, there are two ways by which the 
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parking orbit can be alined with the approach asymptote. Also, there are two ways for 
the parking orbit to shift into alinement with the departure asymptote. In conclusion, 
there are four basic geometries of planetary approach and departure which can be ac- 
commodated by a regressing parking orbit. All four geometries must, of course, be 
scanned to insure finding all available parking orbits. 

Determination of orbital eccentricity . - The orbital eccentricity must now be 

adjusted so that the orbit will rotate the required amount during the planetary stay 

time T . A circular orbit would rotate through the required angles during the 
s 

time T . 
o 

A simple test for the existence of an elliptical orbit that regresses the proper 
amount during the stay time T g is T q < T g . Since circular orbits have the fastest 

regression rates, this inequality means that the eccentricity can be increased to slow 
down the speed of rotation and thereby force T q to equal T . 

The eccentricity that makes T equal to T can be found quite readily by the 

o s 

Newton- Raphson technique. The orbital eccentricity is a real root of a fourth-order 
polynomial found by expressing equation (9) or (10) in terms of the eccentricity. The 
derivation of this polynomial and the solution for orbital eccentricity follow. 

The rate of secular variation of the node and periapsis vectors has been given by 


ft 

s 


- 3nJgR 2 
2a 2 (l - e 2 f 


cos i 


(9) 


S 


-3nJ 2 R 2 
2a 2 (l - e 2 ) 2 


(5 . 2. 
2 Sm 1 


2 ) 


( 10 ) 


The rate ratio ft /o» is determined by the inclination since the terms involving 
s/ s 

eccentricity divide out in the ratio. Thus, the eccentricity can be adjusted independ- 
ently to change the rate without destroying a particular value of the ratio. 

Designating the common factor in both equations by ft f , then 


ft 

r 


-3nJ 2 R 2 

2a 2 (l-e 2 ) 2 


(14) 
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and equations (9) and (10) become 


S2 = f 2 cos i 
s r 


w = « , sin i - 2 
s r\ 


(15) 

(16) 


The mean motion n is 


"“■Jr 3 


(17) 


and the semimajor axis a is 


a = 


1 - e 


(18) 


By substituting equations (17) and (18) for n and a into equation (14) and sim- 
plifying, the following result is obtained: 


n = ( - * VH J,R“r - ' ■ ' ‘ “ I (1 - e) 3/2 (l + e)" 2 


3 . — T _2 -7/2 

r V 2^ J 2 Rr p 


(19) 


3 2-7/2 

For a given radius of periapsis, the quantity - ^ VP- JgR r^ ' is constant. If 


„ 3 / __ „2 -7/2 

c = - 2VP J 2 R r p 


(20) 


then 


K = (1 - e) 3/2 (l + e)" 2 


( 21 ) 
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where 


n 


k = • 


( 22 ) 


Given K, the problem now is to find e. Squaring both sides of equation (21) and 
simplifying, the result is 


K 2 (1 + e) 4 = (1 - e) 3 


(23) 


Expanding equation (23) in a binomial expansion and collecting common powers 
of e results in 


4 (4K 2 + l) 3 (6K 2 - 3 ) 2 (4K 2 + 3) K 2 - 1 ft 

3 s — - e e s e + = 0 


(24) 


K 


K 


K 


K 


Equation (24) will have a real root on the interval 0 < e < 1, provided T < T . Also, 

o s 

equation (24) is of the form 


«/ \ 4 3 2 n 

f(x) = X + CjX + CgX + CgX + C 4 = 0 


(25) 


and the derivative with respect to x of equation (25) is 


f’ (x) = 4x 3 + 3c jX 2 + 2CgX + Cg = 0 


(26) 


The constant coefficients in equations (25) and (26) are 


4K 2 + 1 


_ 4K 2 + 3 
C 3 = tj -2 


C 2 " 


6K 2 - 3 
K 2 


K 


c - k2 - 1 
4 K 2 


> (27) 
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and the constant K is easily evaluated as 


V a s a “n/ T s 

C C cos i C cos i 


(28) 


/ 

where C is the quantity specified in equation (20). 

The parking-orbit eccentricity is found by iteration using equations (25) and (26) 
in 


e 1 = e 
n+1 n 




(29) 


Since e is between 0 and 1, the iteration using equation (29) quickly converges when 
an initial estimate of e = 0. 5 is used. 

A parking-orbit anal ysis program . - The technique discussed has been pro- 
gramed, and the program is essentially a search routine. It accepts two vectors, 

one for the approach trajectory and one for the departure trajectory. The vectors 

define a range of permissable values of inclination for the approach and departure tra- 
jectories. 

The program scans all the possible inclinations and then uses the previously dis- 
cussed interpolation scheme and a Newton-Raphson technique to determine the charac- 
teristics of those parking orbits which match the approach and departure trajectories. 

In general, more than one feasible orbit exists. The program determines all 
parking-orbit configurations .which match the arrival and departure trajectories. The 
output consists of the Keplerian elements of the parking orbit, the approach hyperbola, 
and the departure hyperbola. The orbital parameters of the parking orbit are refer- 
enced to the time of periapsis passage on the approach hyperbola. 


EXAMPLE — 1977 MARS ORBITAL MISSION 


The following example is a preliminary analysis of the parking-orbit phase of a 
near-minimum-energy Mars-orbital mission beginning in 1977. The example illus- 
trates the characteristics of regressing parking orbits which are typical of long- stay- 
time Mars-orbital missions. It also shows how the required characteristics of the 
Mars parking orbit change as the launch date is varied through the Earth-orbital launch 
window. No comparisons are made with alternative methods for parking-orbit aline- 
ment; however, supporting data are included so that evaluations can be made. 
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The interplanetary trajectories and the hyperbolic excess velocities were calcu- 
lated by an independent program. The Earth-orbital launch dates, the trip times, and 
the Mars-orbital stay times were selected so that the heliocentric trajectories would 
be near-Hohmann transfers. These data are shown in figure 3. The Earth-to-Mars 
trip time varied between 360 and 320 days. The Mars-orbital stay time (300 days) and 
the return trip time (320 days) were held constant through the Earth-orbital launch 
window. The total trip time thus varied between 980 and 940 days. The variation of 
the magnitudes is shown in figure 4. 



Figure 3. - Variation of trip times during the Earth-orbital launch window for a 
low-energy Mars-orbital mission beginning in 1977. 


The variation of the AV requirements is shown in figure 5. Earth-orbital 
launch begins from a circular parking orbit at an altitude of 262 n. mi. The parking 
orbit is coplanar with the Earth departure hyperbola. The Mars -orbital insertion 
(MOI) and transearth injection (TEI) impulsive velocities are, of course, dependent on 
the type of parking orbit which is used at Mars. The MOI and TEI impulsive velocities 
are shown for the highest energy regressing parking orbit (orbit 1, table I), which 
shifts into proper alinement for departure. (The orbital geometry code is presented 
in table EL ) The Earth entry velocity (—38 000 fps) is indicated for an entry altitude of 
400 000 feet. 
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Figure 5. - Variation of AV requirements. 


RESULTS 


Ten parking orbits were found to shift into proper alinement for departure. The 
characteristics of these orbits are shown in table I. All the orbits have periapsis alti- 
tudes of 200 n. mi. It should be noted that these orbits represent a set of the highest 
energy parking orbits which shift into proper alinement for departure. There are other 
sets of lower energy orbits which also shift into proper alinement. These sets are 
found by adding integral multiples of 2tt to the required angles of rotation of the 
orbital node Aa 0 and argument of periapsis A cu . In equations (11) and (12), inte- 

P 

gral multiples of 2ir can be added to either Aa^ or Aw^ or both. The resultant 

parking orbits can thus be forced to complete a full revolution or multiple revolutions 
of either the orbital node or the periapsis vector prior to final alinement with the de- 
parture asymptote. 

Table I summarizes the parking orbits available at the beginning, middle, and 
end of the 50-day mission window. The parking orbits are classified into four cate- 
gories corresponding to the four basic geometries of planetary approach and departure. 
If £2^ and £2^ are used to designate the two configurations of the parking orbit ac- 
cording to equation (1) and figure 2, then the four geometries can be designated as 
shown in table II. 


These orbits exhibit several interesting properties. The orbital inclinations fall 
into narrow bands which are widely distributed over the interval of possible inclina- 
tions. Also, the orbital eccentricity is correlated with the orbital inclination. In gen- 
eral, parking orbits with low inclinations have higher eccentricities than more highly 
inclined parking orbits. 


The orbital inclinations are determined in essence by the geometry of planetary 
approach and departure and by the ratio of the required angles of rotation of the orbital 
node and periapsis vectors Aa Q and Aw . The orbital inclination is determined so 


£2 -- p 

that the ratio of the orbital rotation rates 


£2 




u> s is equal to the ratio of the required 


angles of rotation Aa^y'Aw^. The result is that when the required angles of rotation 


are equal and in opposite directions, the rotation rates must likewise be equal and in 
opposite directions. This can occur at inclinations of 46. 4° and 106. 8°. Similarly, 
when the required angles of rotation are equal and in the same direction, the rotation 
rates must be equal and in the same direction. This situation can occur at inclinations 
of 73. 2° and 133. 6°. It would not be surprising, then, if the inclinations of regressing 
orbits fall into narrow regions centered on or somewhere near these equal-rate incli- 
nations. This is the case illustrated in figure 6, which shows the bands of inclination 
of regressing orbits superimposed on a plot of the nodal and periapsis rotation rates. 
For an arbitrary V M geometry, the inclination bands can be shifted to either side of 


the equal-rate inclinations. In cases where the orbit is required to shift through inte- 
gral multiples of 2 tt prior to final alinement, the bands will occur at different charac- 
teristic values of inclination. 
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It was pointed out that parking orbits 
with low inclinations have higher eccentric- 
ities than more highly inclined parking 
orbits. This result can be seen intuitively 
in figure 6. The figure indicates that near- 
equatorial orbits characteristically have 
faster rotation rates than near-polar orbits 
and that their speed of rotation is less 
dependent on eccentricity. Thus, if two & 
parking orbits with identical nodal and per- > 
iapsis rotation rates — one a near-polar 
orbit and the other a near-equatorial 
orbit — are considered, the near- 
equatorial orbit will have the larger eccen- 
tricity. For example, if two orbits have 
inclinations of 46. 4° and 106. 8° (the incli- 
nations for which the rotation rates are 
equal and opposite), the 46. 4° -inclination 
orbit has higher rotation rates than the 
106. 8° -inclination orbit. Thus, the eccen- 
tricity of the 46. 4°-inclination orbit must 
be increased to slow its rates enough to 
match those of the 106. 8 “-inclination orbit. 

Further evidence of the behavior is 
shown in figure 7. Figure 7 shows the 
variation of the orbital eccentricity with 
stay time for a hypothetical, non-varying 



Figure 6. - Secular rotation rates of the 
parking-orbit node and periapsis vec- 
tors due to planetary oblateness. 
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Figure 7. - Variation of orbital eccentricity with stay time for a 
stationary geometry. 



geometry. Three orbits were found to accommodate this geometry. As the stay 

time is shortened, accelerated rotation rates are required to insure final alinement 
with the departure asymptote, and progressively smaller orbital eccentricities are 
required to obtain the faster rates. Figure 8 indicates that a 145° inclined circular 
orbit would rotate through the required angles during a stay time as low as 20 days. 

The curves of figure 8 also show that orbits which are more highly inclined have cor- 
respondingly lower eccentricities. 

Figures 8 and 9 show the characteristics of the highest energy regressing park- 
ing orbit which shifts into proper alinement for departure during the example mission. 
This is the parking orbit for which the MOI AV and TEI AV are given in figure 5, and 
corresponds to orbit 1 of table I. The dashed lines of figure 8 are the inertial traces of 
the orbital groundtrack for the first revolution after MOI and the final revolution prior 
to TEI 300 days later. This is the configuration of the parking orbit as it exists at the 
beginning of the mission window. The solid bell- shaped curve shows the position of 
periapsis at the time of MOI and TEI and at 50-day intervals during the orbital stay 
time. Also, the apparent path of the Sim around Mars is shown at 50-day intervals 
between MOI and TEI. Figure 8 illustrates several interesting characteristics of this 
parking orbit. 



Planetocentric longitude, a, deg 


Figure 8. - Planetocentric motion of periapsis and the subsolar point for a typical 

high-energy parking orbit. 
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Nov. 2, 1977 




Nov. 2, 1977 


Sept. 13, 1977 


Nov. 2, 1977 



3460 


Date of launch from Earth orbit, days after Julian date 2 440 000 


Figure 9. - Required characteristics of a regressing parking orbit and their 
variation during the 1977 mission window. 
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1. The parking-orbit periapsis sweeps through a wide range in latitude, from 
-12° to +18. 6° and then back to -18. 6°. 

2. The periapsis crosses the Mars equator twice — about 40 days after MOI and 
again at 200 days after MOI. This fact would be crucial to minimize plane change re- 
quirements, should an equatorial parking orbit be desired, since the plane change could 
then be performed at the parking-orbit apoapsis. 

3. The periapsis ’’chases" the Sun. Thus, periapsis will remain in sunlight for 
the entire planetary visit. 

Figure 9 shows how the required characteristics of this high-energy orbit vary 
across the Earth-orbital launch window. Basically, the changing characteristics of 
this parking orbit reflect changes in the basic geometry. 


CONCLUDING REMARKS 


Assuming the hyperbolic excess velocities and the periapsis altitude to be non- 
varying, the eccentricity of the parking orbit is the fundamental quantity which governs 
the expense of the Mars-orbital insertion and transearth- injection impulsive velocities. 
The impulsive velocities are greatest for a circular parking orbit and diminish as the 
orbital eccentricity increases. From an impulsive-velocity standpoint, it is therefore 
advantageous to have the parking-orbit eccentricity as large as possible. 

However, orbital eccentricity also governs the speed of rotation of the orbital 
node and periapsis position vectors. Circular orbits experience the fastest regression 
rates. For an orbital mission in which a large plane change must be realized during a 
short stay time, it may be necessary to adopt a near-circular parking orbit. In that 
case, less impulsive velocity may be required to insert into a high- eccentricity ellipse 
and make a parking-orbit plane change. 

The important conclusions regarding the use of this technique therefore rest 
with the trade-off on orbital eccentricity. Obviously, any technique which requires a 
smaller eccentricity will have greater fuel costs. For long planetary stay times (300 to 
500 days), this technique appears to be feasible for orbital eccentricities as large 
as 0. 7 with periapsis altitudes of 200 n. mi. These orbits have inclinations on the 
order of 20°. Higher orbital eccentricities can be obtained by reducing the periapsis 
altitude or the orbital inclination. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, April 16, 1968 
981-30-10-00-72 
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TABLE I. - SUMMARY OF REGRESSING PARKING ORBITS WHICH SHIFT INTO ALINEMENT FOR 


DEPARTURE FOR THREE EARTH-ORBITAL LAUNCH DATES 
(a) September 13, 1977 


. > 1 ; i * i 

Orbit Orbit Inclination, Eccentricity j MOI AV, TEI AV, MOI AV + Total mission 

geometry number deg (a) fps fps TEI AV, fps AV, fps 


1 

1 

18. 64 

0.6974 

3397 

3991 

7 388 

20 041 


2 

68.28 

.5611 

3987 

4581 

8 568 

21 222 


3 

113.36 

.4740 

4379 

4972 

9 351 

22 005 


4 

140. 92 

.6114 

3768 

4361 

8 129 

20 783 

2 

5 

72. 18 

0. 3257 

5071 

5667 

10 738 

23 391 

3 

6 

70. 14 

0.5557 

4011 

4607 

8 618 

21 272 


7 

112. 10 

.4715 

4391 

4985 

9 379 

22 029 


8 

141. 10 

. 6360 

3661 

4255 

7 916 

20 570 

4 

9 

70.95 

0. 3324 

5039 

5633 

10 672 

23 326 


10 

95.04 

. 5105 

4213 

4807 

9 020 

21 674 


a Periapsis altitude = 200 n. mi. 



TABLE I. - SUMMARY OF REGRESSING PARKING ORBITS WHICH SHIFT INTO ALINEMENT FOR 


DEPARTURE FOR THREE EARTH-ORBITAL LAUNCH DATES - Continued 

(b) October 9, 1977 


Orbit 

Orbit 

Inclination, 

Eccentricity 

MOI AV, 

TEI AV, 

MOI AV + 

Total mission 

geometry 

number 

deg 

(a) 

fps 

fps 

TEI AV, fps 

AV, fps 

1 

1 

19. 79 

0. 7018 

3209 

3788 

6 997 

18 823 


2 

68.11 

.5503 

3865 

4443 

8 308 

20 134 


3 

113.30 

.4744 

4207 

4785 

8 992 

20 818 


4 

141.55 

.6125 

3592 

4170 

7 762 

19 588 

2 

5 

59. 56 

0.4518 

4311 

4891 

9 202 

21 027 

3 

6 

70. 13 

0. 5446 

3891 

4469 

8 360 

20 186 


7 

111.86 

.4713 

4221 

4799 

9 020 

20 846 


8 

141.88 

. 6417 

3466 

4044 

7 510 

19 336 

1 4 

9 

58.29 

0. 4480 

4329 

4906 

9 235 

21 061 

1 . 

10 

70. 94 

.3242 

4909 

5487 

10 396 

22 222 


a 


Periapsis altitude = 200 n. mi. 



TABLE I. - SUMMARY OF REGRESSING PARKING ORBITS WHICH SHIFT INTO ALINEMENT FOR 


DEPARTURE FOR THREE EARTH-ORBITAL LAUNCH DATES - Concluded 

(c) November 2, 1977 


Orbit 

geometry 

■l 1 v 

Orbit 

number 

Inclination, 

deg 

Eccentricity 

(a) 

MOI AV, 
fps 

1 

TEI AV, 
fps 

MOI AV + Total mission 
TEI AV, fps AV, fps 

1 

1 

18.46 

0. 7045 

3429 

3646 

7 075 

19 690 


2 

67. 84 

.5513 

4093 

4311 

8 404 

21 018 


3 

113.57 

.4691 

4463 

4681 

9 144 

21 758 


4 

141. 39 

. 6041 

3861 

4079 

7 940 

20 554 


1 2 

5 

71.82 

0. 3245 

5139 

5357 

10 496 

23 110 

3 

6 

69.89 

0.5501 

4101 

4317 

8 418 

21 032 


7 

111.86 

.4655 

4479 

4697 

9 176 

21 790 


8 

141. 72 

. 6384 

3712 

! 

3930 

7 642 

20 256 

4 

9 

58.44 



4765 

9 315 

21 929 



71.24 



5340 

10 466 

23 080 


a Periapsis altitude = 200 n. mi. 






TABLE II. - ORBITAL GEOMETRY CODE 


Orbital geometry 
of regression 

Orbital configuration 
relative to the 
approach asymptote 

1 

°1,A 

2 

°1,A 

3 

°2,A 

4 

°2,A 


Orbital configuration 
relative to the 
departure asymptote 


a 


1) D 


a 


2, D 


a 


1,D 


0 . 


2,D 
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